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ABSTRACT 



The motivation for the search for OvPfi decay is briefly reviewed. It is stressed 
that the exchange of light Majorana neutrinos is not the only possible mechanism 
of the decay. The link between lepton number and lepton flavor violation is 
described and its role in elucidating the Oz//3/3-decay mechanism is discussed. 
Mh| The main topic of the talk is the evaluation of the nuclear matrix elements and 

their uncertainty. Various physics effects that influence the value of the matrix 
elements are described and the results of the two main methods, the quasiparticle 
random phase approximation and the nuclear shell model, are compared. 



> 

G^ , 1. Introduction 

in 

,__! . In the last decade neutrino oscillation experiments have convincingly and tri- 

h^ ! umphantly shown, using both the natural and manmade neutrino sources, that neu- 

O \ trinos have a finite mass and that the lepton fiavor is not a conserved quantity. These 

Q \ results opened the door to what is often called the "Physics Beyond the Standard 

Model". In other words, accommodating these findings into a consistent scenario 
requires generalization of the Standard Model of electroweak interactions that postu- 
lates that neutrinos are massless and that consequently lepton fiavor, and naturally, 
also the total lepton number, are conserved quantities. 

The present results of the oscillation experiments are summarized in Fig. [1] that 
shows the decomposition of the fiavor eigenstates neutrinos z/e, v^j, and v^- into the 
mass eigenstates z/i, z/2 and v^. An upper limit on the masses of all active neutrinos 
~ 2 - 3 eV can be derived from the combination of analysis of the tritium beta-decay 
experiments and the neutrino oscillation experiments. Combining these constraints, 
masses of at least two (out of the total of three active) neutrinos are bracketted by 
10 meV < m^ < 2 - 3 eV. 

Therefore, neutrino masses are six or more orders of magnitude smaller than the 
masses of the other fermions. Moreover, the pattern of masses, i.e. the mass ratios 
of neutrinos, is rather different (even though it remains largely unknown) than the 
pattern of masses of the up- or down-type quarks or charged leptons. All of these 
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Figure 1: Schematic illustration of the decomposition of the neutrino mass eigenstates Ui in terms 
of the flavor eigenstates. The two hierarchies cannot be, at this time, distinguished. The small 
admixture of Ve into v^ is an upper limit, and the mass square of the neutrino vi, the quantity m\, 
remains unknown. 



facts suggest that, perhaps, the origin of the neutrino mass is different than the origin 
(which is still not well understood) of the masses of the other charged fermions. 

The smallness of the neutrino masses can be understood following the finding of 
Weinberg ^ who pointed out almost thirty year ago that there exists only one lowest 
order (dimension 5, suppressed by only one inverse power of the corresponding high 
energy scale A) gauge-invariant operator given the content of the standard model. 
After the spontaneous symmetry breaking when the Higgs acquires vacuum expecta- 
tion value that operator represents the neutrino Majorana mass which violates the 
total lepton number conservation law by two units. 
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where v ~ 250 GeV and C^^-* is expected to be of the order of unity. For sufficiently 
large scale A neutrinos masses are arbitrarily small. 

The most popular explanation of the smallness of neutrino mass is the see-saw 
mechanism, which is also roughly thirty years old ^ . In it , the existence of heavy 
right-handed neutrinos Nji is postulated, and by diagonalizing the corresponding mass 
matrix one arrives at the formula 
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where the Dirac mass m/5 is expected to be a typical charged fermion mass and M^ is 
the Majorana mass of the heavy neutrinos Nj^. Again, the small mass of the standard 
neutrino is related to the large mass of the heavy right-handed partner. Requiring 
that mj, is of the order of 0.1 eV means that M^ (or A) is ~ 10^^"^^ GeV, i.e. near 
the GUT scale. That makes this template scenario particularly attractive. 

Clearly, one cannot reach such high energy scale experimentally. But these sce- 
narios imply that neutrinos are Majorana particles, and consequently that the total 
lepton number should not be conserved. Hence the tests of the lepton number con- 
servation acquires a fundamental importance. 

There are various ways to test whether the total lepton number is conserved or not. 
Examples of the potentially lepton number violating (LNV) processes with important 
limits are 

{Z, A) -^ {Z + 2, A) + 2e-; halflife > 10^^ years 

/i^ + (Z, A) ^ e'^ + {Z — 2, A); exp. branching ratio < 10^"*^^ , 

K~^ — > /i^/i"'"7r~; exp. branching ratio < 3 x 10~^ , 

Ue emission from the Sun; exp. branching ratio < 10~^ . (3) 

However, detailed analysis suggests that the study of the Oi//3/3 decay, the first on the 
list above, is by far the most sensitive test of LNV. In simple terms this is caused by 
the amount of tries one can make. A 100 kg Oi>f3f3 decay source contains ~ 10^^ nuclei 
that can be observed for a long time (several years). This can be contrasted with the 
possibilities of first producing muons or kaons, and then searching for the unusual 
decay channels. The Fermilab accelerators, for example, produce ~ 10^° protons on 
target per year in their beams and thus correspondingly smaller numbers of muons 
or kaons. 

2. Basic considerations 

Double beta decay (/9/3) is a nuclear transition [Z, A) ^ (Z -\- 2, A) in which two 
neutrons bound in a nucleus are simultaneously transformed into two protons plus 
two electrons (and possibly other light neutral particles). This transition is possible 
and potentially observable because nuclei with even Z and N are more bound than 
the odd-odd nuclei with the same A = N + Z. Analogous transition of two protons 
into two neutrons are also, in principle, possible in several nuclei, but phase space 
considerations give preference to the former. 

There are two basic modes of the (3f3 decay. In the two-neutrino mode {2vf3l3) 
there are 2 z/g emitted together with the 2 e~. It is just an ordinary beta decay of two 
bound neutrons occurring simultaneously, since the sequential decays are forbidden by 
the energy conservation law. For this mode, clearly, the lepton number is conserved 
and this mode of decay is allowed in the standard model of electroweak interactions. 
It has been repeatedly observed in a number of cases and proceeds with a typical 




Figure 2: By adding loops involving only standard weak interaction processes the Ov^P decay 
amplitude (the black box) implies the existence of the Majorana neutrino mass. 

half-life of ~ 10^^~^°years. In contrast, in the neutrinoless mode (0z//3/3) only the 2e~ 
are emitted and nothing else. That mode clearly violates the law of lepton number 
conservation and is forbidded in the standard model. Hence, its observation would 
be a signal of a " new physics" . 

One can separate the two modes experimentally by measuring the sum energy 
of the emitted electrons with a good energy resolution, even if the decay rate for 
the Oi//3/3 mode is much smaller than for the 2t//3/3 mode. This is possible since in 
the 2z//3/3 mode the spectrum is continuous, peaked below the Q/2 value of the sum 
kinetic energy, while in the 0z//3/3 mode the spectrum is a 5— function at E1 + E2 = Q, 
smeared only by the experimental energy resolution (the nuclear recoil is always 
negligibly small). 

The existence of the Oz//?/? decay would mean that on the elementary particle level 
a six fermion lepton number violating amplitude transforming two d quarks into two 
u quarks and two electrons is nonvanishing. As was first pointed out by Schechter 
and Valle^ more than twenty five years ago, this fact alone would guarantee that 
neutrinos are massive Majorana fermions (see Fig. [2]). This qualitative statement 
(or theorem), however, does not in general allow us to deduce the magnitude of the 
neutrino mass once the rate of the Oz//3/3 decay have been determined. It is important 
to stress, however, that quite generally an observation of any total lepton number 
violating process, not only of the Oz//3/3 decay, would necessarily imply that neutrinos 
are massive Majorana fermions. 

3. Mechanism of the OuPP decay 

The rather conservative assumption of how the Oz//5/5 decay proceeds is to believe 
that the only possible way the Oz//3/3 decay can occur is through the exchange of a 
virtual light, but massive, Majorana neutrino (the neutrinos z/i, 1/2 and u^ of the first 
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Figure 3: All these symbolic Feynman graphs potentially contribute to the Oi^/3/3-decay amplitude 



section) between the two nucleons undergoing the transition, and that these neutrinos 
interact by the standard left-handed weak currents. 

If we accept this we can relate the 0i//3/3-decay rate to a quantity containing 
information about the absolute neutrino mass. With these caveats that relation can 
be expressed by a well known formula 
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where G^^{Q, Z) is a phase space factor that depends on the transition Q value and 
through the Coulomb effect on the emitted electrons on the nuclear charge, and that 
can be easily and accurately calculated (a complete list of the phase space factors 
^^^(Q,^) and G'^^'^Q^Z) can be found, e.g. in Ref. ^), M^"" is the nuclear matrix 
element that can be evaluated in principle, although with a considerable uncertainty 
and is discussed in detail later, and finally the quantity {mpp) is the effective neutrino 
Majorana mass, representing the important particle physics ingredient of the process. 
In turn, the effective mass {rnpis) is related to the mixing angles Oij (or to the ma- 
trix elements \Ui^i\ of the neutrino mixing matrix) that are determined or constrained 
by the oscillation experiments, to the absolute neutrino masses mi of the mass eigen- 
states Vi and to the totally unknown additional parameters, as fundamental as the 
mixing angles Oij, the so-called Majorana phases a(i). 



{^pp) 



iSJt/. 



|2e*"(*)m,- 



(5) 



Here Uei are the matrix elements of the first row of the neutrino mixing matrix. 

But that is not the only possible mechanism. LNV interactions involving so far 
unobserved much heavier (~ TeV) particles can lead to a comparable Oi//3/3 decay 
rate. Some of the possible mechanisms of the elementary dd —>■ uu + e~e~ transition 



(the "black box" in Fig. [2]) are indicated in Fig. [31 Only the graph in the upper left 
corner would lead to the usual relation between the decay rate and neutrino mass, eq. 
HI Thus, in the absence of additional information about the mechanism responsible 
for the Oz//3/3 decay one could not unambiguously infer the magnitude of {mpj^) from 
the Oi//3/3-decay rate. 

In general Oz//3/? decay can be generated by (i) light mass ive Majorana neutrino 
exchange or (ii) heavy particle exchange (see, e.g. Refs.'^^), resulting from LNV 
dynamics at some scale A above the electroweak one. The relative size of heavy 
{Ah) versus light particle {Al) exchange contributions to the decay amplitude can 
be crudely estimated as follows^: 

where (jripj^) is the effective neutrino Majorana mass, {k"^) ~ (100 MeV)^ is the typical 
light neutrino virtuality, and A is the heavy scale relevant to the LNV dynamics. 
Therefore, Ah/Al ~ 0(1) for {mpp) ~ 0.1 - 0.5 eV and A ~ 1 TeV, and thus the 
LNV dynamics at the TeV scale leads to similar 0i//3/9-decay rate as the exchange of 
light Majorana neutrinos with the effective mass {mpp) ~ 0.1 — 0.5 eV. 

Obviously, the lifetime measurement by itself does not provide the means for de- 
termining the underlying mechanism. The spin-flip and non-flip exchange can be, in 
principle, distinguished by the measurement of the single-electron spectra or polariza- 
tion (see e.g. ^). However, in most cases the mechanism of light Majorana neutrino 
exchange, and of heavy particle exchange, cannot be separated by the observation of 
the emitted electrons. Thus one must look for other phenomenological consequences 
of the different mechanisms. Here I discuss the suggestiorl^ that under natural as- 
sumptions the presence of low scale LNV interactions, and therefore the absence of 
proportionality between {mpp)"^ and the Oz//3/3-decay rate also affects muon lepton 
flavor violating (LFV) processes, and in particular enhances the /i — *> e conversion 
compared to the yU — > 67 decay. 

The discussion is concerned mainly with the branching ratios B^^^,^ = r(yU — >■ 
e7)/r|f) and -B^j^^e = Tconv/rcapt, where /x -^ 67 is normalized to the standard muon 
decay rate V'^^' = (G%ml) / {I927r^) , while /x — *> e conversion is normalized to the muon 
nuclear capture rate F^apt- The main diagnostic tool in our analysis is the ratio 

71 = B^^e/B^^e-f , (7) 

and the relevance of our observation relies on the potential for LFV discovery in the 
forthcoming experiments MEG ' (/i -^ 67) and MECO ^^ (/i ^ e conversionjfj. 



^Even though MECO experiment was recently cancelled, proposals for experiments with similar 
sensitivity exist elsewhere. 



The important quantities are the scales for both LNV and LFV. If they are well 
above the weak scale, then one would not expect to observe any signal in the forth- 
coming LFV experiments, nor would the effects of heavy particle exchange enter Oi>f]f3 
at an appreciable level. In this case, the only origin of a signal in Oz//3/3 at the level 
of prospective experimental sensitivity would be the exchange of a light Majorana 
neutrino, leading to eq.(jl]), and allowing one to extract {mpp) from the decay rate. 

In general, however, the two scales may be distinct, as in SUSY-GUT ' or SUSY 
see-saw ^^ models. In these scenarios, both the Majorana neutrino mass as well as 
LFV effects are generated at the GUT scale. The effects of heavy Majorana neutrino 
exchange in Oi//3/3 are, thus, highly suppressed. In contrast, the effects of GUT- 
scale LFV are transmitted to the TeV-scale by a soft SUSY-breaking sector without 
mass suppression via renormalization group running of the high-scale LFV couplings. 
Consequently, such scenarios could lead to observable effects in the upcoming LFV 
experiments, but with an 0{q) suppression of the branching ratio B^^e relative to 
B^^e-y due to the exchange of a virtual photon in the conversion process rather than 
the emission of a real one. 

The case where the scales of LNV and LFV are both relatively low (~ TeV) is 
more subtle and requires more detailed analysis. This is the scenario which might 
lead to observable signals in LFV searches and at the same time generate ambiguities 
in interpreting a positive signal in Oi>(3f3. This is the case where one needs to develop 
some discriminating criteria. 

Based on the analysis in Ref. '^, we can formulate the main conclusions regarding 
the discriminating power of the ratio 7?.: 

1. Observation of both the LFV muon processes fi —>■ e and fi —>■ e'j with rela- 
tive ratio TZ ~ 10"^ implies, under generic conditions, that Tq^/^ij ~ (m^^)^. 
Hence the relation of the Oz//3/3 lifetime to the absolute neutrino mass scale is 
straightforward. 

2. On the other hand, observation of LFV muon processes with relative ratio 
TZ ^ 10^^ could signal non-trivial LNV dynamics at the TeV scale, whose 
effect on Qvf3(3 has to be analyzed on a case by case basis. Therefore, in this 
scenario no definite conclusion can be drawn based on LFV rates. 

3. Non-observation of LFV in muon processes in forthcoming experiments would 
imply either that the scale of non-trivial LFV and LNV is above a few TeV, 
and thus Tq^pj^ ~ (m/3/3)^, or that any TeV-scale LNV is approximately flavor 
diagonal (this is an important caveat). 

4. Nuclear matrix elements 

Let us assume that the active neutrinos z/i, z/2 and z/3 are indeed massive Majorana 
fermions. If that is so then the neutrinoless (3(3 decay will occur and its rate will be 
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Figure 4: Nuclear matrix elements for the 2i//?/3 decay extracted from the measured half-lives. 



governed by eq.(j4]). Thus, we need to know the value of the nuclear matrix elements 
M^" in order to plan and interpret the experiments. The nuclear transition involved 
consists of changing two neutrons, bound in the ground state (always J'" = 0+) of the 
initial even-even nucleus, into two protons bound in the ground state (again always 
J'^ = O"*") of the final nucleus. (Here we do not consider the case of the excited final 
nuclear states.) 

In the 2v(3(3 decay the two decaying neutrons are uncorrelated. The corresponding 
momentum transfer from the initial neutron to the final proton is small ( q ~ MeV, 
the momentum of the e~ + z/g, qR <^ 1) and thus the long wavelength approximation 
is valid. Since the isospins of the initial and final nuclei are different, only the Gamow- 
Teller operator remains and the corresponding nuclear matrix element is 
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where \m) is the set of all J'^ = 1+ states in the virtual intermediate nucleus. Once 
the lifetime of the 2vl3(3 decay is known, the matrix element M^^ can be easily 
extracted since I/T1/2 = C^^iQ, Z)(Mq^)^, where the phase space factors G'^^{Q, Z) 
can be easily and accurately calculated (seel^). In Fig. H] the M^^ determined in 
this way are depicted. Note the rapid variation of their value when different nuclei 
are involved. 

For the Oz//?/5 decay the situation is quite different. The two initial neutrons, that 
are transformed into the two final protons, are correlated. The virtual massive Ma- 
jorana neutrino that connects them has momentum at least of the order of g > 1/ R 
[R is the nuclear radius, as we will show below the actual momentum transfer is even 
larger) and thus the long wavelength approximation is not valid; all virtual inter- 
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Figure 5: Histogram of older published calculated values of {M'^^)'^ for ^^Ge. The failure of some of 
the calculations to reproduce the known 2zy/3/3-decay lifetime is indicated. 



mediate states can in principle contribute. On the other hand, the typical nuclear 
excitation are less than the energy of the virtual neutrino, hence the "closure ap- 
proximation" is valid, and we usually do not need to worry about the energies of 
intermediate states. 

In nuclear structure theory one begins with the mean field approximation in which 
the nucleons are bound in a potential, but independent. That is, however, a poor 
approximation, and the "residual interaction" need to be taken into account, typically 
with some truncation. There are two basic and complementary methods of evaluating 
M^'^, the quasiparticle random phase approximation (QRPA and its various general- 
izations) and the nuclear shell model (NSM). These two methods differ fundamentally 
in the way the indicated truncation is implemented. In QRPA one selects a wide in- 
terval of single-particle orbits but only a class of configurations (particle-hole and its 
iterations) of the nucleons are taken into account. On the other hand, in NSM only a 
relatively narrow interval of single-particle orbits is chosen (one oscillator shell or less) 
but all (or almost all) configurations of the valence nucleons residing on those orbits 
are included. In both methods an effective interaction is used, based on the known 
nucleon-nucleon force, but modified slightly using selected nuclear data as guidance. 
Since these two methods are so different, it is important to test them against each 
other. 

There are many evaluations of the matrix elements M^'^ in the published literature 
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(for the latest review see^"^). However, the resulting matrix elements often do not 



agree with each other and it is difficult, based on the pubhshed material, to decide 
who is right and who is wrong, and what is the theoretical uncertainty in M^'^. That 
was stressed in the paper by Bahcall et al. few years ago '^^ where a histogram of 20 
calculated values of {M'^^)'^ for ^^Ge was plotted, with the implication that the width 
of that histogram is a measure of uncertainty. That is clearly not a valid conclusion 
as one could see in Fig. [5] where the failure of the outliers to reproduce the known 
2i//5/3-decay lifetime is indicated. 




Figure 6: Contributions of different angular momenta J^ of the two participating neutrons to the 
Gamow-TeUer part of M^'^ in *^Se (upper panel) and ^'^'^Te (lower panel). The results of the 
nuclear shell model (NSM) (dark histogram) and QRPA treatments (lighter histogram) are com- 
pared. Both calculations use the same single-particle spaces: {f5/2iP3/2,Pi/2i 99/2) for *^Se and 
(97 /2, d5/2, d3/2, Si/2, hii/2) for ^■^''Te. In the QRPA calculation the experimental 2i'f3(3-deca.y rate 
was used to adjust the strength of the particle-particle interaction. 

So why are the calculated values of different authors different? In order to under- 
stand the difficulties in evaluating M^'^ we plot in Fig. [6] the contribution of different 
angular momenta J^ of the two transformed neutrons. There are two opposing ten- 
dencies in Fig. [6l The large positive contribution (essentialy the same in QRPA and 
NSM) is associated with the so-called pairing interaction of neutrons with neutrons 
and protons with protons. As the result of that interaction the nuclear ground state is 
mainly composed of Cooper-like pairs of neutrons and protons coupled to J^ = 0. The 
transformation of one neutron Cooper pair into one Cooper proton pair is responsible 
for the J7 = piece in Fig. [61 



However, the nuclear hamiltonian contains, in addition, important neutron-proton 
interaction. That interaction, primarily, causes presence in the nuclear ground state 
of "broken pairs", i.e. pairs of neutrons or protons coupled to ^7 7^ 0. Their effect, 
as seen in Fig. [6l is to reduce drastically the magnitude of M^'^ . In treating these 
terms, the agreement between QRPA and NSM is only semi-quantitative. Since the 
pieces related to the "pairing" and "broken pairs" contribution ale almost of the same 
magnitude but of opposite signs, an error in one of these two competing tendencies is 
enhanced in the final M^". The competition, illustrated in Fig. [6l is the main reason 
behind the spread of the calculated M°^. Many authors use different, and sometimes 
inconsistent, treatment of the neutron-proton interaction. 

In our QRPA (and RQRPA, renormalized QRPA) calculations l ^^ l ^'^ l ^^^ we adjust 
the neutron-proton particle-particle interaction, responsible for the "broken pairs" 
contribution, using the known 2z//3/3-decay lifetimes, respectively the corresponding 
M^^ matrix elements. We based this procedure on the fact that the Gamow- Teller 
strength, the contribution of the 1+ virtual intermediate states that is fully responsible 
for the M"^" , is the quantity most sensitive to the corresponding parameter, usually 
denoted as Qpp. The nominal value, corresponding to the G-matrix based on the 
realistic nucleon-nucleon force, is gpp = 1. The renormalized gpp has values between 
0.8 and 1.2. 

In the NSM there is no analog to adjustment of the Qpp parameter. The hamil- 
tonian (primarily its so-called monopole part) is adjusted so that a set of nuclear 
spectroscopic data is optimally reproduced. That is a complicated and time consum- 
ing procedure, but as a result a number of nuclear properties is well described. There 
is no attempt to reproduce specifically the 2 i//j/3-deca y lifetime, but the agreement 
with experiment is, in most cases, acceptable l -i-^ l ^u i /i ii 

To see some additional reasons why different authors obtain in their calculations 
different nuclear matrix elements we need to analyze the dependence of the M^'^ on 
the distance r between the pair of initial neutrons (and, naturally, the pair of final 
protons) that are transformed in the decay process. That analysis reveals, at the 
same time, the various physics ingredients that must be included in the calculations 
so that realistic values of the M^" can be obtained. 

The corresponding Ou^P decay operator contains, besides the spin and isospin 
operators, the "neutrino potentials", i.e. the Fourier transform of the neutrino prop- 
agator. The simplest and most important of these potentials has the form 

H(r) = -<^(ujr) , (9) 

r 

where R is the nuclear radius introduced here to make the potential, and the resulting 
M'^'", dimensionless (the 1/R^ in the phase space factor compensates for this), r is 
the distance between the transformed neutrons (or protons) and ^{ur) is a rather 
slowly varying function of its argument. 




Figure 7: The dependence on r of M^" for '^^Ge, ^^^Mo and ^^°Te. The upper panel shows the 
full matrix element, and the lower panel shows separately 'pairing' [J = of the two participating 
neutrons) and 'broken pair' [J ^ 0) contributions. 



From the form of the potential H{r) one would, naively, expect that the charac- 
teristic value of r is the typical distance between the nucleons in a nucleus, namely 
that f r^ R. However, that is not true as was demonstrated first in Ref. ^^ and 
illustrated in Fig. [71 One can see there that the competition between the "pairing" 
and "broken pairs" pieces essentially removes all effects of r > 2 — 3 fm. Only the 
relatively short distances contribute significantly. The same result was obtained in 
the NSM l^^. (We have also shown in ^^ that an analogous result is obtained in 
an exactly solvable, semirealistic model. There we also showed that this behaviour 
is restricted to an interval of the parameter gpp that contains the realistic value near 
unity.) 

Once the r dependence displayed in Fig. [7] is accepted, several new physics effects 
come to mind. One of them is the short-range nucleon-nucleon repulsion known 
from scattering experiments. Two nucleon strongly repel each other at distances 
r < 0.5 — 1.0 fm, i.e. the distances very relevant to evaluation of the M^'^ . The 
nuclear wave functions used in QRPA and NSM, products of the mean field single- 
nucleon wave function, do not take into account the influence of this repulsion that is 
irrelevant in most standard nuclear structure theory applications. The standard way 
to include the effect is to modify the radial dependence of the Qujdp operator so that 
the effect of short distances (small values of r) is reduced. 
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Figure 8: The r dependence of M^'^ in ^^Ge. The four curves show the effects of different treatments 
of short-range correlations. The resulting M"'' values are 5.32 when the effect is ignored, 5. Of when 
the UCOM transformation is applied and 4.14 when the treatment based on the Fermi hypernetted 
chain and 3.98 when the phenomenological Jastrow function is used. (See the text for details.) 



Traditionally, this was done by multiplying; the neutrino potential by the square 
of a Jastrow-like function first derived inl^ and in a more modern form inl^. That 
phenomenological procedure reduces the magnitude of M^'^ by 20-25% as illustrated 
in Fig. [HI Recently, another procedure, based on the Unitary Correlation Operator 
Method (UCOM) has been proposed 1^. That procedure, still applied not fully con- 
sistently, reduces the M^'^ much less, only by about 5%'^. It is prudent to include 
these two possibilities as extremes and the corresponding range as systematic error. 
Once a consistent procedure is developed, consisting of deriving an effective OuPP 
decay operator that includes (probably perturbatively) the effect of the high momen- 
tum (or short range) that component of the systematic error could be substantially 
reduced. 

Another effect that needs to be taken into account is the nucleon finite size. That 
is included, usually, by introducing the dipole form of the nucleon form factor 



/y,. 



i + qyui^^ 



(10) 



where the cut-off parameters My^A have values (deduced in the reactions of free neu- 
trinos with free or quasifree nucleons) ~ 1 GeV. This corresponds to the nucleon 
size of ~ 0.5 — 1.0 fm. Note that in our case we are dealing with neutrinos far off 
mass shell, and bound nucleons, hence it is not obvious that the above form fac- 
tors are applicable. It turns out, however, that once the short range correlations are 
properly included (by either of the procedures discussed above) the M^'^ becomes 
essentially independent of the adopted values when Mv,a > 1 GeV. However, in the 
past various authors neglected the effect of short range correlations, and in that case 
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Figure 9: The full ranges of M '^^ with the two alternative treatments of the short range correlations. 
For comparison the results of a recent Large Scale Shell Model evaluation of M '^'" that used the 
Jastrow- type treatment of short range correlations are also shown (triangles). 



a proper inclusion of nucleon form factor (or their neglect) again causes variations in 
the calculated M'^^ values. 

Yet another correction that various authors neglected must be included in a correct 
treatment. Since r < 2-3 fm is the relevant distances, the corresponding momentum 
transfer 1/r is of the order of ~200 MeV, much larger than in the ordinary /3 de- 
cay. Hence the induced nucleon currents, in particular the pseudo scalar (since the 
neutrinos are far off mass shell) give noticeable contributions l -i-o l -^Qt . 

We have, therefore, identified the various physics effects that ought to be included 
in a realistic evaluation of M^" values. The spread of the calculated values, noted 
by Bahcall et al. ^^ can be often attributed to the fact that various authors either 
neglect some of them, or include them inconsistently. 



5. Discussion and conclusion 

Even though we were able to explain, or eliminate, a substantial part of the 
spread of the calculated values of the nuclear matrix elements, sizeable systematic 
unce rtainty remains. That uncertainty, within QRPA and RQRPA, as discussed in 
RefsJ ^ ^ , is primarily related to the difference between these two procedures, to the 
size of the single-particle space included, whether the so-called quenching of the axial 
current coupling constant qa is included or not, and to the systematic error in the 
treatment of short range correlations^^. In Fig. [9] the full ranges of the resulting 
matrix elements M^'^ is indicated. The indicated error bars are highly correlated; 
e.g., if true values are near the lower end in one nucleus, they would be near the lower 
ends in all indicated nuclei. 



The figure also shows the most recent NSM results '^^. Those results, obtained 
with Jastrow type short range correlation corrections, are noticeably lower than the 
QRPA values. That difference is particularly acute in the lighter nuclei ^^Ge and ^^Se. 
While the QRPA and NSM agree on many aspects of the problem, in particular on 
the role of the competition between "pairing" and "broken pairs" contributions and 
on the r dependence of the matrix elements, the disagreement in the actual values 
remains to be explained. 

When one compares the 2z/ and Oz/ matrix elements (Figs. H] and [9]) the feature to 
notice is the fast variation in M'^'^ when going from one nucleus to another while M^'^ 
change only rather smoothly, in both QRPA and NSM. This is presumably related 
to the high momentum transfer (or short range) involved in Oi>f3(3. That property of 
the M^^ matrix elements makes the comparison of results obtained in different nuclei 
easier and more reliable. 

While substantial progress has been achieved, we are still somewhat far from 
being able to evaluate the Oz//?/3 nuclear matrix elements confidently and accurately. 
Perhaps at the next Venice meeting we will be able to report that we reached that 
goal. 
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